Large spin splitting of GaN electronic states induced by Gd doping 
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We present a detailed study of the magnetic-field and temperature-dependent polarization of the near-band- 
gap photoluminescence in Gd-doped GaN layers. Our study reveals an extraordinarily strong influence of Gd 
doping on the electronic states in the GaN matrix. We observe that the spin splitting of the valence band reverses 
its sign for Gd concentrations as low as 1.6 x 10 16 cm~ 3 . This remarkable result can be understood only in 
terms of a long range induction of magnetic moments in the surrounding GaN matrix by the Gd ions. 
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The current interest in semiconductor spintronics provides 
a strong motivation for fundamental studies of diluted mag- 
netic semiconductors (DMS) OlEEIl- Apart from the advan- 
tage of making all-semiconducting spintronic devices accessi- 
ble, one particular benefit of DMS is the possibility to control 
their electronic and magnetic properties by the same means 
jj]. Among the different II- VI and III-V semiconductor mate- 
rials, GaN-based DMS have the potential to fulfill the require- 
ment of high-temperature ferromagnetism for future spintronic 
devices |5]. Room-temperature ferromagnetism in transition- 
metal (TM) doped DMS, however, is often found to result from 
the formation of precipitates rather than from a homogeneous 
alloy 1 6]. Rare-earth (RE) elements could be an interesting 
alternative, since they have partially filled / orbitals, which 
carry magnetic moments and may contribute to magnetic cou- 
pling like in the case of TM with partially filled d orbitals. 
Gd is the only RE element which has both partially filled 4/ 
and 5d orbitals. Together, these orbitals can take part in a new 
coupling mechanism proceeding via intra-ion A-f-5d exchange 
interaction followed by inter-ion 5d-5d coupling mediated by 
charge carriers 0]. Teraguchi et al. j|| observed a Curie tem- 
perature larger than 400 K in the Gao.94Gdo.06N alloy. Very 
recently, we have observed ferromagnetism above room tem- 
perature as well as an extraordinarily large magnetic moment 
of Gd in weakly Gd-doped GaN layers |9]. The average value 
of the magnetic moment per Gd ion is found to be as high as 
4000 flB as compared to its atomic moment of 8 }1b. This find- 
ing suggests that the GaN matrix strongly contributes to the 
observed macroscopic magnetization. 

Optical spectroscopy has been demonstrated to be a pow- 
erful experimental tool to study the interplay between elec- 
tronic and magnetic properties in DMS fiol [Tli O . In 
this work, we used magneto-photoluminescence (MPL) spec- 
troscopy to investigate the influence of Gd doping on the elec- 
tronic properties of the GaN matrix. The energy and polariza- 
tion of near-band-gap optical emission are expected to be di- 
rectly influenced by spin splittings in the electronic band struc- 
ture. We studied, in particular, the circular polarization of the 
MPL emission due to recombination of donor-bound excitons 
(D°,X), which typically dominates the low-temperature PL in 
GaN. Our study reveals a strong Gd-induced spin splitting of 
the valence band, which has an opposite sign as compared to 
the Zeeman splitting in these samples. 

The 600-nm-thick GaN:Gd layers were grown with Gd con- 



centrations (Nai) of 1.6 x 10 16 (sample A) and 6 x 10 16 cm -3 
(sample B) directly on 6H-SiC(0001) at a substrate tempera- 
ture of 8 10 °C using reactive molecular-beam epitaxy (RMBE) 
121 . An undoped GaN reference sample was grown under 
identical conditions. The concentration of Gd was determined 
by secondary-ion-mass spectrometry (SIMS) and found to be 
constant over the entire depth. All samples were subject to 
an extensive investigation by high-resolution x-ray diffraction 
over a wide angular range. No reflections related to a sec- 
ondary phase were detected. Furthermore, sample B was in- 
vestigated by transmission electron microscopy. The layer did 
not contain any clusters or precipitates. Electrically, all sam- 
ples, including the reference sample, are found to be highly re- 
sistive. The magnetic properties of a series of Gd-doped GaN 
layers, including samples A and B, have been systematically 
investigated by superconducting-quantum-interference-device 
(SQUID) magnetometry. Colossal magnetic moments of 3000 
and 1000 pLg per Gd ion were found in samples A and B, re- 
spectively. Both samples are also found to exhibit ferromag- 
netism above room temperature. For MPL experiments, the 
samples were excited by the linearly polarized emission of a 
He-Cd ion laser at a wavelength of 325 nm. The laser power 
density focused on the sample was varied between 200 and 
250 Wcm~ 2 . Since the ferromagnetism in the GaN:Gd layers 
is not induced by free carriers (investigated layers are insulat- 
ing) 10], we do not expect a significant change of the mag- 
netic characteristics due to the small density of photocreated 
carriers. The PL light was dispersed in a single-path spectro- 
graph and detected by a charge-coupled device (CCD) array. 
The experiments were carried out in the temperature range 5- 
100 K in a continuous He-flow cryostat and in magnetic fields 
up to 12 T using the backscattering Faraday geometry. In all 
measurements, the magnetic field was oriented parallel to the 
c-axis of the samples. The magnetic field-induced circular po- 
larization of the PL light was analyzed by passing it through 
a photoelastic modulator (PEM), a quarter-wave plate, and a 
linear polarizer. 

In the absence of an external magnetic field, the PL spec- 
tra of all samples are characteristic for undoped epitaxial GaN 
layers in that they are dominated by the (D° 7 X) transition at 
a photon energy of 3.458 eV. The lower energy of this transi- 
tion, when compared to homoepitaxial GaN, is consistent with 
the tensile in-plane strain in these layers of 0.15% lfl3ll . In all 
cases, the donor responsible for this transition is oxygen with 
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FIG. 1 : Circularly polarized PL spectra for (a) the reference sample 
and (b) sample B measured at 7 K under a magnetic field of B = 10 T 
in the Faraday configuration (B || c). 



a concentration of about 10 18 cm~ 3 as measured by SIMS. 
Since these donors are distributed homogeneously over the en- 
tire GaN matrix, the (D° 7 X) emission can be utilized as a probe 
for the properties of the electronic band structure. The forma- 
tion of Gd-O complexes can be neglected since there is only 
one Gd ion available for hundreds of oxygen donors at our 
low Gd-doping concentrations. Figure shows the PL spec- 
tra for an undoped GaN reference sample and sample B under 
an external magnetic field of B = 10 T. The observed (D°,X) 
emission is polarized in both samples, which is evident from 
the difference in intensities of the two circularly polarized a 
(full squares) and a + (open squares) components. Most im- 
portantly, the polarization for sample B has the opposite sign 
of the one in the reference sample. This finding is remarkable 
since it implies that a relatively small amount of Gd is able to 
reverse the sequence of Zeeman splitted states in a large vol- 
ume fraction of the GaN matrix. 

The (D°,X) complex consists of three spin particles, two 
electrons (one each in the donor and exciton) and one valence 
band hole (in the exciton). In the ground state, two electrons 
are in antiparallel spin orientation (singlet state), while the spin 
of the valence band hole is uncompensated. Consequently, the 
circular polarization of the PL emission due to (D°,X) recom- 
bination is determined only by the population of spin-up and 
spin-down states for holes. Note that, based on similar argu- 
ments, acceptor bound excitons (A°,X) could be used to study 
the influence of the Gd doping on electron states. The uncom- 



pensated hole spin has a non-zero projection only along the c- 
axis because the ground state of (D°,X) in GaN is associated 
with the A exciton. An external magnetic field along the c-axis 
thus splits the (D°,X) ground state into two Zeeman sublevels 
with up and down projection of the hole spin. Here we assume 
that the Zeeman splitting (due to an external magnetic field) 
is much smaller than the singlet-triplet energy splitting of the 
electron pair in the complex. Hence, the (D°,X) emission is 
expected to be composed of two types of circularly polarized 
light (a~ and (7 + ). The relative population of the two Zee- 
man sublevels determines the circular polarization degree of 
the emitted light, which is defined as 
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where l a± denotes the intensity of the a ± polarized light, kg 
the Boltzmann constant and T the temperature. To and z s are 
the lifetime of the (D°,X) transition and the spin-relaxation 
time of the holes, respectively. To account for the variation 
of the spin-relaxation time with the magnetic field observed in 
Ref.Ojl we assume 
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where a and / are constants and Bo = 12 T the maximum mag- 
netic field. Note that the spin relaxation time for electrons and 
holes can be expected to be similar because of the weak spin- 
orbit interaction in GaN. Indeed, for semiconductor quantum 
wells and dots, having non-degenerate valence bands like in 
wurtzite GaN, long spin relaxation times between 1 and 20 ns 
have been reported 1 15U16I1 . AE(B) = ghl^-BB is the energy sep- 
aration between the Zeeman sublevels of holes at an external 
magnetic field B. The g-factors for electrons g e (= 1.95 (3) 
and A-valence band holes g/,( = 2.25 1 18]) are quite close to 
each other in GaN. Thus, the splitting between the a and <T + 
polarized (D°,X) emission is expected to be w 0.1 meV at B 
= 10 T, which is too small to be resolved. However, the same 
magnetic field can produce a noticeable polarization of the PL 
emission as apparent from Fig. The sign reversal of the 
polarization in sample B with respect to the reference sample 
(Fig. [0 implies that AE(B) has changed its sign. For a mag- 
netic semiconductor, AE(B) = ghl^BB + AEi nt , where AEi nt is 
the additional splitting generated by the internal field produced 
by the magnetic atoms which can be regarded as a renormal- 
ization of the g-factor. In our case, AEi nt is evidently negative 
and stronger than the Zeeman term, which is possible only 
when the exci tonic ground state is subject to a strong influ- 
ence of the Gd ions. Since the Gd-doped samples also do not 
reveal any energy splitting between the a and <7 + polarized 
(D°,X) emission, we conclude that the conduction and valence 
bands are influenced in the way by the presence of Gd ions (see 
arguments above). 

Figure |2 displays the (D°,X) polarization p as a function 
of magnetic field for sample A (circles), B (squares) and the 
reference sample (triangles). Clearly, the magnetic field de- 
pendence of the polarization for the Gd-doped samples are 
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FIG. 2: Circular polarization p of the (D ,X) emission as a function 
of the external magnetic field B measured at 7 K for the reference 
sample (triangles) as well as samples A (circles) and B (squares). 
Solid lines are the theoretical fit obtained from our model discussed 
in the text. 

different from that of the reference sample. While the refer- 
ence sample exhibits, as expected, a nearly linear magnetic- 
field dependence, a tendency toward saturation is observed in 
particular for the Gd-doped sample B. It is also evident from 
Fig. |2 that the relative change of the polarization with respect 
to the reference sample increases with the concentration of Gd. 
All these findings confirm that a large fraction of excitons is 
influenced by the presence of Gd. 

A sign reversal and a qualitatively similar magnetic-field de- 
pendence of the polarization of the PL emission due to the 
recombination of neutral Mn-acceptor-bound excitons (A°,X) 
has been observed in (Ga,Mn)As which is a III-V DMS. The 
observed behavior of the PL polarization has been explained 
in terms of an antiferromagnetic coupling between the valence 
band hole and the Mn cZ-orbital |19]. In this case, however, 
the excitons are bound to the magnetic dopant itself and are, 
consequently, easily influenced by the Mn of-orbitals. On the 
contrary, the excitons in our Gd-doped GaN layers are bound 
to O donors, which are homogeneously distributed within the 
matrix. The mean Gd-Gd separation in sample B is 25 nm re- 
sulting in an average distance between a (D°,X) site and a Gd 
ion as large as 12 nm 1 20] . None of the existing theories can 
explain a magnetic coupling across such a large distance 12 ill . 

The extraordinary situation encountered here becomes ob- 
vious when we attempt to analyze our data in the conventional 
framework valid for DMS such as (Ga,Mn)As, i. e. under 
the assumption of a direct interaction between Gd ions and 
(D°,X). It is then expected that the additional splitting A£i nt 
follows the external-field dependence of the macroscopic mag- 
netization. SQUID measurements reveal a ferromagnetic be- 
havior in these samples with a saturation of the magnetiza- 
tion above B —2 T in both samples A and B. For the sake of 
simplicity, let us consider sample B in the following example. 
A fit of Eq. Q to p obtained in the magnetization-saturation 
range (cf. Fig. [5} yields A£; nt = —2.5 meV. However, AZ?i nt 
= NoPx/icd, where fled = 8 jJ-s is the bare atomic moment of 
Gd, No is the number of Ga or N atoms per unit volume, x = 
Ned/No, and Nofi is a constant which defines the strength of 



FIG. 3: Circular polarization of the (D ,X) emission as a function 
of the external magnetic field B measured at different temperatures 
for sample B. Solid lines are theoretical fits obtained from our model 
(discussed in the text). 



TABLE I: Parameters obtained from the fits for different GaN:Gd 
samples used in this study: Gd concentration Nqj (10 16 cm -3 ), sam- 
ple temperature T (K), additional splitting A£; nt (meV) generated by 
the Gd ions, parameters a and y for the magnetic-field dependence of 
the hole spin relaxation time, volume fraction y ocuupied by the Gd- 
spheres of influence, and corresponding radius r (nm) of the spheres 
(cf. Eq.gJ. 
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the exchange coupling between the Gd ions and the valence 
band holes. The corresponding value of Afoj3 ~ —200 eV is 
two orders of magnitude larger than the ones obtained in II- VI 
and III-V DMS and thus much too large to be physically mean- 
ingful. A reasonable value of Nofi is only obtained when the 
effective magnetic moment per Gd ion is two orders of magni- 
tude larger than its atomic moment of 8 jUg, which is consistent 
with the colossal magnetic moment of 1000 \Ib obtained from 
the magnetization measurements on sample B. Repeating this 
analyis for sample A, we obtain an even larger value for A^j3, 
in agreement with the fact that the measured magnetic moment 
per Gd ion is as high as 3000 Hb in this sample. In agreement 
with our previous conclusion,]^, we thus have to consider that 
the Gd ions induce a magnetic moment in a large number of 
matrix atoms (Ga and/or N) in order to understand the efficient 
spin-polarization of a large fraction of excitons by the Gd dop- 
ing. 

To account for this polarization of the matrix, we associate 
a sphere of influence with each of the randomly positioned Gd 
ions. Since the radius r of these spheres is presumably much 
larger than the lattice spacing in GaN, we further assume that 
the spheres are randomly arranged in a three-dimensional con- 
tinuum (continuum percolation). The polarization p of the 
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(D ,X) emission can now be expressed as 

P = (1 -^Pout+yPin (3) 

where y is the volume fraction of the regions occupied by the 
spheres of influence. Within the framework of continuum per- 
colation, y can be expressed as 

y = 1 - e- vN ™ (4) 

where v = 47Tr 3 /3 is the volume of the sphere of influence, 
in and out are the contributions to the polarization by excitons 
bound to donors inside and outside the regions spanned by the 
spheres, respectively. Both p; n and p out can be expressed by 
Eq. Q, whereas AE(B) is given by: 

AEout(fi) = ghllBB 

AE in (B) = gkHsB + AEint. (5) 

The values of a and y (cf. Eq. are expected to be dif- 
ferent inside and outside the regions occupied by the spheres 
of influence. z s (B) for regions not occupied by a sphere of in- 
fluence is required to be identical to the one in undoped GaN, 
which can be obtained from a fit to the data for the reference 
sample by Eq. Q. 

Next, we use this model to fit the experimental field depen- 
dence of the polarization at 7 K (Fig.g). The data which are 
obtained in the saturation range of the magnetization {B > 2 T) 
|U, are only taken in to account in this fitting. A£; nt could thus 
be considered as a constant in this range. Only a, y and y are 
used as fit parameters. Note that the hysteresis at B <2 T can- 
not be resolved by our magneto-PL setup, since the coercive 
field (0.01 T) as well as the remanent magnetization are too 
small 1 9]. The central result of this approach is that satisfac- 
tory fits are obtained only with finite values for y, particularly 
so for sample B, for which y has to be very close to unity. A 
volume fraction of y = 1 means that the entire matrix is un- 
der the influence of the Gd ions, i. e., the radius of the sphere 
of influence is close to the average Gd-Gd spacing (25 nm) 
in this sample. Consequently, the fitting procedure was first 
applied to the reference sample and to sample B, representing 
the cases of y = and y=l, respectively. The obtained values 
of a and y were subsequently taken into account to calculate 
the polarization curve of sample A using only A£; nt and y as 
fit parameters. This fit yields r = 17.2 nm, which is close to 
the value obtained for sample B demonstrating the consistency 
of this approach. The relation between the volume fraction y 
and the magnetic moment per Gd ion is explained in detail in 
Ref.H 

A further test of our model is given by the temperature 
dependence obtained for sample B. As shown in Fig. |3j our 
model is able to reproduce the data with reasonable values for 
all fit parameters (Tab. [j}. Furthermore, the temperature de- 
pendence of A£j nt is in reasonable agreement with the decrease 
in saturation magnetization observed in the SQUID measure- 
ments. Finally, we point out that our findings are consistent 
with the size of the Gd sphere (r « 30 nm) estimated in Ref. 9 
to account for the colossal magnetic moment per Gd ion ob- 
served in these samples, which provides further quantitative 



support in favor of a long-range interaction of Gd with the 
GaN matrix. 

Our study reveals an extraordinarily large effect of Gd on 
the electronic states in GaN. It is observed that even a Gd con- 
centration as low as 1.6 x 10 16 cm~ 3 can significantly influ- 
ence the spin splitting of the conduction and valence bands. 
This unprecedented effect can be explained only in terms of a 
long-range induction of magnetic moments in the GaN matrix 
by the Gd ions. This conclusion is consistent with the colos- 
sal magnetic moment per Gd ion and ferromagnetism with 
Tc > 300 K observed in these samples. Our results are, in par- 
ticular, interesting for spintronic applications, since GaN:Gd 
may be easily doped with donors (acceptors) with a concentra- 
tion exceeding that of Gd to generate spin-polarized electrons 
(holes) in the conduction band (valence band). Gd-doped GaN 
with its Curie temperature above room temperature may thus 
be a very attractive candidate as a spin-injector material for 
all-semiconducting devices 12211 . 
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